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We present a dynamical theory of the Auger decay in one-dimensional (ID) metals described 
by the Tomonaga-Luttinger model. An analytic expression of the Auger current is derived in the 
framework of the 1-step approach, where the finite lifetime of the initial core-hole and the core- 
valence interaction are taken into account. This allows to capture typical dynamical features like 
the shake-down effect, in which the Auger spectrum shows a non-vanishing weight above the 2-step 
high-energy threshold. The obtained results give also a hint to understand the sizable suppression 
of Auger spectral weight closed to the Fermi energy recently observed in carbon nanotubes with 
respect to graphite. 



INTRODUCTION 

The dynamics of the core-valence-valence (CVV) 
Auger transitions in strongly correlated solids has been 
extensively studied during the last three decades 
However, despite the great interest devoted to this prob- 
lem, several aspects are still not well understood. The 
theoretical calculation of the Auger spectrum of corre- 
lated solids is a challenging task because, beside the 
intrinsic difficulty to deal with a many-body interact- 
ing system, the creation of the core-hole and the Auger 
process itself are in principle coherent events, involv- 
ing virtual Auger transitions and incomplete relaxation 
phenomena [2]. The complete formulation of the theory 
describing the Auger decays has been provided in 1980 
by Gunnarsson and Schonhammer (GS)[3]- In the frame- 
work of the so-called 1-step approach, they derived a gen- 
eral formula for the Auger current by treating the decay 
of the initial core- hole to all orders. Unfortunately such 
formulation cannot be cast in terms of Green's functions 
and is of hard implementation for practical purposes. A 
significant progress can be done within the 2-step approx- 
imation, where the photoemission and the Auger decay 
are considered as independent events. In this framework, 
CiniQ and SawatzkyQ (CS) proposed a simple and ele- 
gant theory able to provide a quantitative understanding 
of the experimental Auger spectra of transition metals 
with (almost) closed valence bands. An advantage of 
such theory is that it also provides a practical scheme to 
estimate the value of the screened interaction from the 
experimental spectra @, 0]- This is particularly useful 
to support LDA-I-J7 calculations In the case of open- 
band systems the CS approach breaks down and no reli- 
able theory is currently available. Very recently Seibold 
and collaborators presented a theory of the dynamical 
two-particle response function in the 2D Hubbard model 
based on the time-dependent Gutzwiller approximation. 
Although important effects are not treated there (e.g., 
the finite lifetime of the core-hole and the interaction of 
the core- hole with the valence electrons), the theory pro- 
vides a novel tool to attack the calculation of the Auger 



spectrum in correlated open-band systems. 

In this paper we develop a dynamical theory of the 
CVV Auger transitions in an ideal ID metal. In the 
CVV Auger decay two holes are left in the valence band 
after the X-ray photoemission of a deep core-hole. Here 
assume that the valence electrons form the so-called Lut- 
tinger Liquid (LL) , described by the Tomonaga-Luttinger 
model. We also allow for the interaction between the 
core-hole and the valence electrons and introduce a term 
responsible for the Auger transition, which destroys the 
core-hole and creates the Auger electron together with 
the two valence holes (and viceversa). The correspond- 
ing Auger current is calculated analitically by using the 
bosonization and equations of motion methods. 

The paper is organized as follows. In Section H we de- 
scribe in detail the model Hamiltonian used in the present 
work. In Section III we derive a closed analytic expres- 
sion for the Auger current in framework of the 1-step 
approach. In Section IV we discuss the relevant features 
emerging form the obtained formula. In particular we 
show that the theory exposed here is able to capture some 
striking features recently observed in the Auger spectra 
of carbon nanotubes. Finally, a brief summary and the 
the main conclusions are drawn in Section V. 



THE MODEL 

The LL is the prototype of interacting electrons con- 
fined in one spatial dimension, characterized by strik- 
ing phenomena such as the so-called spin-charge sepa- 
ration and the power-law dependence of observables in 
proximity of the Fermi energy [13, U, 12, 13, [13, [l5| . 
In -ffLutt the electrons have a linear dispersion relation 



around positive (Right) and negative (Left) Fermi points 
and the electron-electron (e-e) interactions act only be- 
tween Right/Left electron densities. The model is exactly 
solvable by means of the bosonization technique, which 
allows to write the electron Hamiltonian in terms of bo- 
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where cr =t,i is the spin index, [ba{q),bl^,{q')] — 
Sa;a-'Sq,q' , vp is the Fermi velocity, is the interaction 
parameter between Right-Right (positive q) and Left- 
Left (negative q) electron densities, while (72 is the in- 
teraction parameter between Left-Right densities. 

The key point of the bosonization is that it is possible 
to express the fermion fields in terms of boson fields. Here 
it is useful to introduce a chirality index v = R, L to 
distinguish between Right and Left electron modes. For 
instance the v = R fermion field is given by: 
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where rjR,^ is an anticommuting Klein factor and 



g>0 ' 



27r 



1/2 



aq/2 



bliqy 



(3) 



where is the total number of Right-electrons, 

[Vo,Rj ^r] ~ i ^-iid L is the length of the system, a is a 
short-distance cutoff that must be introduced in order to 



have converging integrals [l7j. In principle the bosoniza- 
tion provides exact results in the limit a — > 0, however for 
practical purposes it is useful to take a non-zero (small) a 
which introduces a finite effective bandwidth 7 = vpja 
in the system. By doing this we have to bare in mind 
that such procedure gives an accurate physical descrip- 
tion only in the low-energy part of the spectrum [l8l|. 

The cou pling of valence electrons to the core-hole is 
given byiUllli^ 



- «c) , (4) 



where A(g) is the core- valence coupling constant, L is the 
volume of the system, ST' is the annihilation (creation) 
operator of the core- electron, whose occupancy and en- 
ergy are ric = c^Cc and £c respectively. In the following 
we will take A(g) = A. 

The term responsible for the Auger decay is more con- 
veniently expressed in the fermionic representation and 



Ha = 44 A + h.c. , A = Vi;^- (0)^x (0) , (5) 
where Cp'' destroys (creates) the Auger electron and 

Ipa- = ^a,R + Ipa.L ■ (6) 

V is the so-called Auger matrix element, which here is 
taken as a constant. Here we are assuming for simplic- 
ity that the CVV decay leaves the two final holes in the 
origin of the system in a singlet configuration. This re- 
fiects the local nature of the Auger process; however such 
assumption is not essential and could be relaxed. 

As long as the interactions do not depend on spin, 
^^Lutt can be diagonalized by introducing charge and 
spin boson operators: fec|l(g) — [6^^^ (q) ± 6|^^ ((7)]/V2 and 
performing a Bogoliubov transformation in the charge 
sector: bc{q) — cosh(pbc{q) + sinhipbl{—q), b\{q) — 
sinh (p6c(— <?) + cosh(p6j((j) with tanh2(p = {g2/'K)/{vF + 
gi/T^) and renormalized velocity v — [{vp + g4/7r)^ — 
(sz/tt)^]^^^. In the next Section we use the bosonization 
scheme sketched above compute the Auger spectrum of 
a ID metal described within the Luttinger liquid theory. 



CALCULATION OF THE AUGER SPECTRUM 

The 1-step formulation of the the Auger processes has 
been provided by GS[2][2l| who showed that the Auger 
current is given by the following correlator [2^ : 



7ra 



2 poo poo 



dt dt' e^-'^'-'"^ J{t,t') (7) 
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where the factor 7ra^/2 is chosen in order to have a nor- 
malized spectrum and 

f{t,t') 

= (g|cte^(«["l+'nt'^^^te^i/[i](t-0^cte-'^(^M-'f)*c,|.g).(8) 

In the above expression \g) is the ground state before the 
X-ray photoemission, H[0, 1] is the Hamiltonian of the 
system H = Hi,^tt+£c{^ — 'nc)+H\ with Uc = 0, 1 respec- 
tively, and r is an effective optical potential describing 
virtual Auger transitions and relaxation processes [2]. In 
order to proceed we make the following approximation: 

f{t,t') « (g|e'^["l*'Ate*^W(*-*')^e-*^["l*|.g) 



X e-*^'e(t-t')g-r(W) 

= C(i,t')xe-*^=(*-*')e-^(*+*') 



(9) 



where the second line of the above equation is the core- 
hole Green's function with lifetime l/F which is a c- 
number. \g) denotes the ground state of -ffLutt (whose el- 
ementary excitations are created by b^ and bl), describing 
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the valence band in the initial state, and H[0] = H[0]—ec- 
As noticed by GS0, the strength of the effective optical 
potential is proportional to the square of the Auger ma- 
trix element, and hence we can replace V"^ by F. C{t,t') 
can be calculated exactly by using the bosonization for- 
mulas in Eqs. (|2l3p and the equations of motion method. 
After some algebra one gets a compact expression by in- 
troducing new variables t — t — t' and T — {t + t')/2: 



r 

2^ 
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(10) 

where we have defined g = 2 (cosh ip + sinh ip) and and 
I = (cosh 99-l-sinh (^)^. The complex functions h(T,T) and 
k(T,T) are reported in Appendix A. Finally the Auger 
current reads 
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where we refer the kinetic energy of the Auger electron 
u! with respect to the core- level energy Ec- Egs. fTU)) and 
pip constitute the main finding of the present work. In 
the next Section we discuss the most relevant features 
emerging from Egs. ljlOp and (jlip . which give a hint to 
understand the physics of the Auger transitions in ID 
systems. 



DISCUSSION 

We first observe that despite the LL nature of the va- 
lence electrons, the correlator C(r, T) does not obey a 
power-law, which is spoiled by the interaction A between 
the valence electrons and the core-hole. It is also inter- 
esting to study the relationship of our solution with the 
2-step approach. This is done in the limit F ^ 0. As 
discussed by GS, if such limit exists, one should recover 
the well-known 2-step solution since the Auger transition 
happens after the complete relaxation of the initial state. 
Such limit is carried out by observing that 
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with z ^ h,k. We note that for any finite A the limit 
on the r.h.s. does not exist because for large T we have 
/i(t,T) ~ fc(r,r) - {X/vfivTlii{vT/a). This is a re- 
markable result, showing that the 2-step approach is not 
justified if the valence band is described by the LL. On 
the other hand if we set A = the 1-step and 2-step solu- 
tions do coincide because the 2-step spectrum is obtained 
from the 2-particle Green's function describing the va- 
lence electrons (holes) in the ground state of 7?Lutt • The 
2-step approach is often employed in typical Auger calcu- 
lations and therefore it is instructive to compare it with 
our 1-step solution. The 2-step Auger current is readily 



FIG. 1: Auger current j{uj — Ec) calculated numerically from 
Ea. (|ll|l for different values of the core-hole lifetime and core- 
valence interaction. Here we have taken g4 = g2 = 4, a = 0.1, 
uf = 1, ffl = 1, and A = 4 except in black curve, where A = 0. 
ij is expressed in units of 7. The vertical dashed line denotes 
the 2-step high-energy threshold Ec, set equal to zero in the 
present figure. 



obtained by setting A = and results 



j2-stcp('i^ - £c) 



dT- 



{—iTV + a)9 (—iTV + ay{—iTVp + a 



which recovers the characteristic power law suppression 
at cj « 0. The comparison between the Auger spectra 
calculated with j{uj — Ec) and j2-stcp(w — Ec) is shown in 
FigU As discussed above, it is seen that j does not ap- 
proach j2-stcp for small F (compared to 7). In particular 
we note that for F ^ the center-of-gravity Eg of j (blue 
and violet curves) is shifted towards lower kinetic ener- 
gies with respect to the center of gravity of j2-stcp (black 
curve) with a logarithmic dependence Eg oc A^ ln(Fa/ti). 

Another interesting case is obtained in the limit F — > 
00, that is for very short core-hole lifetime. In this case 
2Fe~^'""^ produces a Dirac delta in the T-integration: 



hm 2F 



(14) 



with z — h,k. In this limit the Auger transition occurs 
when the initial state is still excited, since e"*^^*'!*!^) is 
not an eigenstate of H[l]. As a consequence the exci- 
tations created on emission of the initial core-electron 
transfer their energy to the Auger electron, which then 
has a kinetic energy exceeding the high-energy thresh- 
old Ec in the 2-step model [see Eq. (fT3|) and FigU] (black 
curve)]. This phenomenon, known as shake-down, is a 
typical example of qualitative departures from the pre- 
dictions of the 2-step model. 
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Finally, we show that our theory provides an expla- 
nation of the sizable suppression of CVV Auger spec- 
tral weight closed to the Fermi energy observed in car- 
bon nanotubes with respect to graphite [13 [H. This IS 
a striking trend, since the structure of the one-particle 
density of states (IPDOS) of the two carbon structures 
would predict the opposite behavior. In fact, while in 
metallic nanotubes the IPDOS at the Fermi energy is 
finite due to the ID linear dispersion, in graphite it is 
vanishing, due to the 2D conical dispersion. Therefore 
we expect that correlation effects have to be invoked in 
order to revert the one-particle scenario. 

Metallic carbon nanotubes are believed to be 
rather good (although approximate) realizations of 
LLji^ljlallil] since in normal conditions the main cor- 
relation effects come from the long-range part of the 
Coulomb repulsion. In nanotubes with radius R, the 
back-scattering interactions with large momentum trans- 
fer suffer a 1/i? suppression [ill [l^. Therefore we believe 
that typical metallic (10,10) nanotubes are well described 
by the present theory. Concerning graphite, we use the 
CS approach, which is known to give the Auger spectrum 
in excellent agreement with experimental one 2J| 30| . In 
order to employ the CS approach we must compute the 
2-particle valence Green's function within the bare ladder 
approximation Q This is accomplished starting from 
the non-interacting valence IPDOS 



Ah 



(15) 



which is obtained by imposing the 2D linear spectrum 
eikx,ky) = vp{k'^ + fcp^/^ and the momentum cutoff 
1/a. We note en passing that vanishes linearly at 
w = 0, as it should. The corresponding non-interacting 2- 
particle Green's function Gq^ is obtained by self-folding 
and by Hilbert transforming: 



G^^(^) ^ (l/67)[2-(c./7) + 2(c./7)^ 
- (a;/7)^e-/^r(0,c./7)]. 



(16) 



where r(x, y) is the incomplete gamma function. Thus 
the Auger spectrum of graphite according to CS theory 
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where U is the short-range screened repulsion felt by the 
two valence holes in the final state. It is worth to recall 
that our model is suitable to represent the tt electrons of 
nanotubes and graphite, which are the ones involved in 
proximity of the Fermi level (placed at a; = 0). Therefore 
only the low-energy portion of the experimental Auger 
spectra can be addressed within the present framework, 
while the high-energy spectral region, corresponding to 
deep Up and as states, cannot be described here. A com- 
plete analysis including the missing Up and Gs states can 
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FIG. 2; Comparison between the calculated Auger current for 
graphite jcs (^ ^ £c) (dashed black line) and (10,10) SWNTs 
jiuj — Ec) according to Eq. (|11|) (solid lines) for different values 
of the core-valence interaction, uj is expressed in units of 7 
and Ec is set equal to zero. 



be found in Ref.[2j]. However, in that paper the sup- 
pression of the Auger spectrum of nanotubes closed to 
the Fermi energy was reproduced by including some phe- 
nomenological form factors which have been fitted with 
the experimental data. Conversely in the present work 
the problem is treated starting form a fully microscopic 
theory with no adjustable parameter. 

In order to compare with the experiment, we use the 
following realistic values for graphite and (10,10) metal- 
lic nanotubes: vf ~ lO^m/s, a such that 7 = vf/ol ~ 
lOeV, F « 0.2eV[33LXi.e^ F/7 = 0.02), =_a2.~ 



2e2 \n{L/R)l 




and leaving the adimensional ratio X/^/aj as free pa- 
rameter. In Figl5] we see that the inclusion of e-e cor- 
relations in carbon nanotubes according to the LL the- 
ory within the 2-step approach (black solid line) is not 
enough to reproduce the suppression of j vs j^g* closed 
to w ~ 0. On the other hand the full 1-step theory with 
finite A/ ^/a ^ 7 provides results in qualitative agreement 
with the experimental trend [s^. This is a quite reason- 
able finding, considering that the core-valence repulsion 
is larger than (but of the same order of) the valence- 
valence repulsion. 



SUMMARY AND CONCLUSIONS 

Traditional photoemission and inverse photoemission 
experiments which probe one-particle dynamical re- 
sponses provide a well-established tool for the under- 
standing of strongly correlated ID systems. A great 
amount of theoretical work devoted to this problem has 
been published in the past, enlightening the role of the LL 
concept to explain several features 3, 11 , 13, 13, Id , lB| ■ 



Surprisingly the study of the ID Auger transitions, 
which are related to the two-particle dynamics, has been 
only poorly addressed. On the other hand the Auger 
spectroscopy is a powerful experimental technique which 
permits the characterization of the correlations in solids 
and therefore is of crucial importance in the study of 
strongly correlated systems. 



In the present work we have developed a dynamical 
theory of the Auger processes in ID metals described 
within the LL theory. Our theory includes the finite core- 
hole lifetime, the valence- valence and the core- valence in- 
teractions as well. A typical 1-step feature is observed in 
the limit of small core-hole lifetime, in which the valence 
electrons cannot relax before the Auger transition, and 
the shake-down phenomenon occurs. Remarkably it is 
shown that the 2-step approximation is not valid for any 
finite core-valence interaction, which also spoils the low- 
energy power-law behavior typically expected in the LL. 
Only for vanishing core-valence interaction the power- 
law is recovered. Finally we have shown that our 1-step 
theory is able to reproduce the low-energy suppression 
of Auger spectral weight observed in carbon nanotubes 
with respect to graphite. 
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discussions. 
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with A = A\/2(cosh ip — sinh p). 



APPENDIX A: THE FUNCTIONS /i(r,T) AND 
k{r,T) 



The final expression for C{t,t') in Eq. (fTO|) has been 
obtained by employing the bosonization formulas in 
Eas. (|2l3[ ) and the equations of motion method. In or- 
der to perform the sum over q we used that q = 2Trn/L 
and took the large-L limit. When doing this, it is useful 
to set L = aN where N is the number of sites of the ID 
system and a is the lattice constant, and send — > oo. 
The functions /i(r, T) and fc(r, T) obtained in this way 
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